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Abstract: We analyze the existence, stability, and propagation of 
dissipative discrete localized modes in one- and two-dimensional nonlinear 
lattices composed of weakly coupled split-ring resonators (SRRs) excited 
by an external electromagnetic field. We employ the near- field interaction 
approach for describing quasi- static electric and magnetic interaction 
between the resonators, and demonstrate the crucial importance of the 
electric coupling, which can completely reverse the sign of the overall 
interaction between the resonators. We derive the effective nonlinear model 
and analyze the properties of nonlinear localized modes excited in one- 
and two-dimensional lattices. In particular, we study nonlinear magnetic 
domain walls (the so-called switching waves) separating two different 
states of nonlinear magnetization, and reveal the bistable dependence of the 
domain wall velocity on the external field. Then, we study two-dimensional 
localized modes in nonlinear lattices of SRRs and demonstrate that larger 
domains may experience modulational instability and splitting. 
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1. Introduction 

Important building blocks of electromagnetic metamaterials |1| are the split-ring resonators 
(SRRs) or other types of subwavelength resonant elements which are arranged in one-, two-, or 
three-dimensional lattices. In general, the response of a metamaterial is not simply given by a 
sum of the responses of individual resonators, but it depends also on the near-field interaction 
between the resonators within the system ||2l[3l . A standard theoretical approach for analyzing 
the properties of metamaterials is based on the effective medium approximation when the struc- 
ture is treated as a homogeneous medium characterized by effective macroscopic parameters. 
This approximation is justified when the characteristic wavelength of the electromagnetic field 
is much larger than the period of the micro structured medium. However, for smaller wavelength 
or in the cases when the internal resonant modes are important, the metamaterials demonstrate 
their discrete nature and strong nonlocal effects [4], so they should be described as lattices of 
resonant elements by employing the techniques usually used in the analysis of photonic crystals 
and waveguide arrays |i5J . 

The discrete physics of composite metamaterials can be studied using a novel type of coupled 
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Fig. 1. Electric and magnetic interaction coefficients (solid and dashed curves, respectively) 
as functions of the spacing between the ring centers for the normalized frequency = 
—0.2, when the rings are offset in (a) x direction, (b) y direction, and (c) z -direction. Insets 
show schematically corresponding ring positions. 
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Fig. 2. Examples of (a,b) arrays and (c,d) two-dimensional lattices of weakly coupled non- 
linear split-ring resonators creating the simplest nonlinear magnetic metamaterials. 



nonlinear equations which describe the so-called nonlinear magnetoinductive waves [6] when 
the magnetic response of a metamaterial may become bistable. In this paper, we extend the 
earlier analysis of Ref . O and consider a lattice of SRRs analyzing various types of nonlinear 
dissipative discrete localized modes. Although the nonlinear discrete equations were originally 
derived in Ref. [6| and further analyzed in Refs. iTTl-fTTIl, here we introduce an important gen- 
eralization of those equations which takes into account both electric and magnetic near-field 
coupling between the neighboring sites, as was recently described by a deeper analysis of res- 
onator interaction ||2|[T2l . 

We demonstrate that in the nonlinear regime the magnetic response of an array of split-ring 
resonators may become bistable, so that a nonlinear metamaterial can support the propagation 
of domain walls (also called switching waves) separating the regions of different values of the 
magnetization, and study the motion of such domain walls under the action of an external field, 
revealing the hysteresis in the dependence of the velocity on the applied field. We also analyze 
the existence and stability of nonlinear localized modes in two-dimensional magnetic lattices 
of SRRs including the decay of magnetic domains via modulational instability. 

The paper is organized as follows. Section 2 is devoted to the derivation of the new nonlinear 
model that takes into account both electric and magnetic interaction between the split-ring 
resonators. This new model is then applied in Sec. 3 to the study of nonlinear localized modes in 
dissipative one- and two-dimensional lattices composed of nonlinear SRRs. More specifically, 
we confirm the existence of magnetic domain walls in arrays of SRRs when both the interaction 
terms are taken into account (Sec. 3), and also study the nonlinear dynamics of dissipative 
localized modes in two-dimensional lattices when modulational instability may break up the 
broader modes into smaller modes (Sec. 4). Finally, Sec. 5 concludes the paper. 

2. Near-field interaction between split-ring resonators 

The equation describing the current In^q^m in a split-ring resonator with indices (n,q,m), corre- 
sponding to the axes (x,y,z), has the form O 

^ = (mj^ +Pjej) , (1) 

where L,7?,C are inductance, resistance and capacitance of the resonator, Qn^q^m is the charge 
in the resonator, j describes indices of the nearest neighboring resonators, M and P are electric 
and magnetic interaction coefficients, respectively, (on^q^m is an external electromotive force. 

Following 1 6], we assume that the capacitance of the resonator is nonlinear, i.e. it depends 
on the voltage U across the gap of the resonator: Cnl = Q + AQv^ (|^w,^,m)P), where the 
nonlinear correction to the capacitance is small. Using the slowly varying approximation (61 
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Fig. 3. (a) Dependence of the magnetization on the index nfoYL = 0.024. Sohd - standing, 
and dashed - moving domain wall, corresponding to the top and bottom branches on the 
panel (c) in this Figure, (b) Dependence of the velocity of the domain wall as a function of 
the external excitation E. (c) shows close up of the bistable region in panel (b), with arrows 
indicating directions of the jumps. All results are obtained for Q. = —0.2, = —0.02. 



we obtain the following equation for the normalized current in the split-ring resonator: 

where the dimensionless variables T = cOot, cOq = I/LCq, = {co — C0o)/c0o, and ^n,q,m = 
^n,q,ml where J^c = cOqCqUc is the characteristic nonlinear current, Uc = Ec x dg is the 
characteristic nonlinear voltage, 7 = R/LcOq is the damping coefficient, and 

is the normalized electromotive force, ^n,q,m is the amplitude of the harmonic current In^q^m- 
The effective interaction coefficients K-^ = k^h^ / ^ — ^,e^ / where normalized electric 
and magnetic interaction coefficients are k^ h = M/L and k^ e = P/LcOq. These interaction 
coefficients crucially depend on the mutual position of the resonators, and the distributions of 
currents and charges in conductors . 

To be more specific, in further calculations we assume that the resonators have the single ring 
geometry identical to that used in Ref. O, namely: rings have average radius ro=2.25 mm, track 
width of 0.5 mm, metal thickness of 0.03 mm, gap width of 1 mm. For such resonators, the in- 
teraction coefficients were calculated using our approach presented in Ref. [2J, where we have 
also taken into account the effect of retardation [13]. The calculated coefficients are presented 
in Fig. ([T]). We verified that the resonant frequencies of a pair of interacting rings described by 




Fig. 4. Examples of a nonlinear localized mode in a two-dimensional lattice of SRRs. a = 
1,^2 = —0.2, 7= 0.024, E = 0.058. (a) soliton with minimum possible size (b) soliton 
with size 5x1. 



the interaction coefficients match the resonances found from direct numerical simulations us- 
ing commercial software CST Microwave Studio. This makes us confident that the coefficients 
can quantitatively describe near-field interaction of the resonators. Importantly, the effective 
interaction coefficient in Eq.[2]depends on the difference between electric and magnetic co- 
efficients, which for experimentally realistic parameters may become comparable. In contrast 
to the previous works, which assumed only approximate magnetic interaction between the res- 
onators, when the interaction of in-plane resonators is identical in x- and y-directions, in our 
case we see from the Figs.[T](a,b) that the interaction is strongly different in x- and y-direction. 
Remarkably, the corresponding effective interaction coefficients K can have different signs. 

3. Discrete dissipative localized modes 

3.1. Domain walls in SRR arrays 

First, we study switching wave propagation in one-dimensional arrays of split-ring resonators, 
see Figs. [3l extending our previous analysis lfT4l . We select the geometry shown in Fig. [2a) 
with the spacing between the resonators of c/ro = 2.4. From Fig. [T] we find that the corre- 
sponding effective interaction coefficient = —0.02. In the bistable regime, when two stable 
uniform solutions exist in the chain, we observe the possibility of formation of switching waves, 
which represent a transition from one uniform distribution to another with a change in resonator 
number n. Such waves are initially excited by non-uniform distribution of the external field, and 
then they are supported by a uniform external excitation. The profiles of such waves are close 
to a step function, with typical domain wall structure shown in Fig.Oa). 

Because of the discreteness of the system, the velocities and profiles of moving switching 
waves vary quasi-periodically with time (in contrast to a switching wave with a stationary 
profile in continuous media). For this reason, the velocity of the switching wave v is obtained 
by averaging over sufficient time. The velocity sign is defined positive if the motion of the 
switching wave leads to the expansion of the region occupied by the upper branch of bistability ; 
otherwise, the velocity is negative. Figure [3tb,c) shows the velocity of the switching waves as 
a function of the external excitation. We find that for the selected parameters there exist a 
narrow range of external excitations L for which we observe bistable behavior, see Fig. He). 
Interestingly, one of the branches in bistable regime corresponds to the stationary domain wall, 
while another branch - to the one moving in negative direction. 

3.2. Dissipative solitons in two-dimensional lattices 

Now we study the two-dimensional arrays of the resonators, and we select the lattice orientation 
shown in Fig.^c). We choose the normalized frequency of the excitation of ^1 = —0.2, and the 
spacings between resonators a/ro = 2.4, c/ro = 1. From the results presented in Fig. [T] we find 




Fig. 5. (a) Initial and (b) final state of the two-dimensional magnetic domain in a SRR lattice 
in the case of anisotropic modulational instability. The domain decays in two identical 
modes of a smaller size. Parameters are: a = I, Q. = —0.2, 7= 0.032, E = 0.063. 



the effective interaction coefficients for this regime as /T^ = 0.06, = —0.02, and use these 
values in our further simulations of Eqs. [21 Our two-dimensional structure supports a family 
of nonlinear localized solutions - discrete dissipative solitons. In bistable regime, when the in- 
teraction between resonators is sufficiently weak, such solutions of Eqs. [2] can be found using 
perturbation theory, with the smallness parameter defined by the interaction coefficients. How- 
ever for realistic interaction constants, shown in Fig.[TJ the existence and stability of dissipative 
solitons can be found only numerically. Fig. |4] shows two possible stable soliton solutions, con- 
sisting of just one strongly excited resonator (Fig.Ua)) or five excited resonators (Fig.|4tb)). 
For such narrow solitons, as well as for narrow solitons in continuous media flS], stability of 
the upper branches of the bistability curves is not generally required for the stability of the soli- 
tons, since the modulational instability is calculated for the homogeneously excited structure. 
However, for wider excitations, the modulational instability starts to manifest itself, with the 
larger structures showing instability and splitting into several smaller solitons. For example, if 
we initially excite the structure with wide area corresponding to the upper branch of the bista- 
bility curve (see Fig. [3a)), then the excitation becomes unstable in z-direction, and we observe 
formation of two stable narrow solitons, which are shown in Fig. [3b). Demonstrated decay of 
the wide soliton occurs within 100 dimensionless time units T. 

4. Conclusions 

We have analyzed nonlinear modes in one-dimensional arrays and two-dimensional lattices 
composed of weakly coupled split-ring resonators with both electric and magnetic coupling. 
We have demonstrated the existence of domain walls connecting two stable states of the meta- 
material excited by an external electromagnetic field, and studied its dissipative dynamics char- 
acterized by the bistable velocity. We have also studied two-dimensional nonlinear magnetic 
domains in nonlinear lattices of SRRs and demonstrated that larger domains may experience 
modulational instability. Our results may be useful for the analysis of other types of nonlinear 
metamaterials supporting subwavelength discrete localized modes, especially those operating 
in the optical regime, such as nanoscale periodic structures consisting of metal and nonlinear 
dielectric slabs where discrete solitons have been analyzed recently |[T6ll . 
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